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Geometry and velocity structure of HD 44179's bipolar jet 
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ABSTRACT 

In this paper we analyse a set of 33 optical spectra, which were acquired with 
the ARCES echelle spectrograph (R = 38,000) on the 3.5-ni telescope at the Apache 
Point Observatory. We examine the Hq profile in each of these observations in order 
to determine the geometry and velocity structure of the previously discovered bipolar 
jet, which originates from the secondary star of HD 44179 located at the centre of 
the Red Rectangle nebula. Using a 3D geometric model we are able to determine the 
orbital coverage during which the jet occults the primary star. During the occultation, 
part of the Ha line profile appears in absorption. The velocity structure of the jet was 
determined by modelling the absorption line profile using the Sobolev approximation 
for each orbital phase during which we have observations. The results indicate the 
presence of a wide angle jet, likely responsible for observed biconical structure of the 
outer nebula. Furthermore, we were able to determine a likely velocity structure and 
rule out several others. We find that the jet is comprised of low-density, high-velocity, 
central region and a higher-density, lower-velocity, conical shell. 
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1 INTRODUCTION 

The Red Rectangle (RR) is an X-shaped, biconical, pro- 
toplanetary nebula associated with the star HD 44179. Its 
moniker is derived fr om its rectangular appearance in red 
photographic plates l| Cohen et al.lll975l ). while its appear- 
ance in the blue plates is approximately circular. The ap- 
pearance of the RR in the red is dominat ed by its high 
intensity of extended red emission (ERE) (|Schmidt et al.1 
ll980l : IWitt fc Boroson|[i990 ). powered by the central source 
(HD 44179), which is a single line spectroscopic binary 
van Winckel et al]|l995l '). The binary consists of a post- 



AGB primary star (0.8 Mq) with an effective temperature 
of roughly 8000 K, and a near solar mass main-sequence 
companion that has an accretion disk and jet (see Witt et 
al. 2009, and references therein). 

The X-shaped structure corresponds to the walls of a 
biconical low-density cavity within the spherical AGB out- 
flow. The walls of the cavity appear bright where the line of 
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sight is tangential to the cone. The walls of the cavity are 
directly illuminat ed by the cen tral source , which g ives rise to 
the strong ERE (|Schmidt fc W itt I99l1: ICoheii et al . 20o3; 
IViih et al.i[200^ ). The far U V photons (energ ies > 10.5 eV) 
necessary to excite the ERE (|Witt et al.ll200g ) are produced 
in th e hot accretion d isk (T ~ 17,000 K) of the secondary 
star JWitt et al.ll2009l ). 



The RR's seemingly unique appearance can be at- 
tributed to its orientation on the sky with respect to our 
view from Earth. The central binary is obscured by a near- 
edge-on optically-thick circumbina ry disk with an i nclina - 
tion angle of 86° as measured by Buiarrabal et al.l (|2005l ) 
or 85° as measured by iThomas et al.1 l|201ll ) Therefore, the 
central source can only be observed via scattering over the 
top and bottom of the circumbinary disk. This indirect line 
of sight causes us to view the binary at an effective inclina- 
tion ang le of 35° (|Waelkens et al.lll996l ). The effective incli- 
nation is an average angle, as the orbital motion produces a 
variation of ± 0.8°. Since the inferred properties of the sys- 
tem are dependent upon knowing the correct angle for our 
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line of sight to the primary star, all such derived properties 
are somewhat uncertain. 

Significant research has been conducted on the morphol- 
ogy of the extended nebula. The most likely agent for pro- 
ducing the morphology is a high speed bipolar jet launched 
from the accretion disk around the secondary star. The bipo- 
lar jet proceeds to carve a low density cavity in the more 
slowly expanding spherical cloud of material ejected during 
the primary star's AGB stage. Th is scena r io for shapin g such 
bipolar nebulae was proposed bv lMorrid (|l98ll . ll987l ). Such 
a high s peed jet has been observed in the RR via Ha ab- 
sorption (|Witt et al.ll2009l ). Studies of the Ha spectra reveal 
that the velocity of th e jet material rang es between a min- 
imum of 150 k m s~^ (iKoning et aL l bOllI ) 
of 560 km s~ 
Recent 



Witt et aLlbOoi) 



and a maximum 



work by Velazquez et al.l (|201ll ) and 



iKoning et al.l (I2OIII ) have used the binary model to 
explain other feature s of the morphology of the RR. 
IVelazguez et al.l (|201lh posit a narrow precessing jet that 
originates from the secondary star with the goal of produc- 
ing the lad der-like ru ngs that appear in the outer nebula. 
Converselv. IKoning et a l. (2011) employed a wide angle jet 
to evacuate the biconical cavity. Using a scattered light 
model (with the biconical cavity devoid of dust) they show 
that the ladder rungs could be ring-like structures that 
are observed in other bipolar nebulae, which in the case of 
the RR, when pr ojected onto the plane of the sky, appear 
as ladder rungs. ISokeil l|2005h demonstrates that special 
physics is not required to shape the RR, and suggests 
intermittent jets give rise to the ladder like structures in 
the RR. 

The goal of this paper is explore the geometry and kine- 
matics of the jet using an analysis of the Ha spectra as a 
function of orbital phase. Is the jet narrow or wide? Specif- 
ically, we demonstrate that the opening angle of the jet 
is roughly equal to the geometric opening angle of the X- 
shaped structure of the RR. Therefore, our data are more 
consistent with the model of iKonirig et al.l (|201ll ). 

In this paper, we first introduce the framework in which 
the data will be discussed. The orbit and the physical model 
for the jet are described in § [2] We then present the obser- 
vations in § O and the measurements in § [J] Calculations 
based on a numerical model and a discussion are presented 
in §[5] We present our conclusions in §[6] 



2 MODEL 

A 3D geometric model of the RR system is necessary for 
interpreting our measurements (to be discussed in § 2} . The 
principal components of this system are: the primary star, 
the secondary star and its accretion disk, the bipolar jet, the 
circumbinary disk, and the spherical outflow. We approxi- 
mate the jet as a bi-cone centr ed on the sec ondary star, and 
include the orbital solution of I Thomas et al. (2011) for cal- 
culating the positions of the various objects as a function of 
orbital phase. 



2.1 Geometry 

The inclination angle, i, of the RR is 85°. The distance 
to the RR remains uncertain, but we adopt the value 



of 710 pc derived by iMen'shchikov et all (|2002l ') for our 
calculations. The outer radius of the circumbinary disc 
is 1850 au (|Buiarrabal et al"] |2005| ). while the radius of 
the central cav ity in the circumbinary d isk is estimated 
to be 14 au (|Men'shchikov et all 120021 '). The thickness 
of the circumbi nary disk is determined t o be approxi- 
mate l y 100 au jMen'shchikov et al.l l2002l: iRoddier et all 



I1995I : lOsterbart et all Il997l : iBuiarrabal et all bOOSl ). The 

opening angle (angle from cavity wall to cavity wall) of 
the biconical cavity is 40° just above the circumbinary 
disk, while at larg er distances it expands from 62° to 80° 
(|Cohen et al.ll2004 ). For modeling purposes we use 70° (cor- 
respond ing to a 35° hal f open ing angle), as was used in the 
work of IWaelkens et al.l (|l996l ). 

Figure [1] illustrates the innermost region of our model 
for the RR. The outer radius of the accretion disk (~0.3 au) 
is set to be the distance of the LI Lagrange point from the 
secondary star when the binary components are at perias- 
tron. The size of the primary star (radius ~ 0.21 au) and 
the orbital separation of the binary components are shown 
to-scale. Notice that the direct line of sight (direction to 
the observer) is blocked by the optically thick circumbinary 
disk. Consequently, the primary star is viewed indirectly (at 
an effective inclination angle, ies = 35°) via scattering over 
the top and bottom of the circumbinary disk. As a result, 
the indirect line of sight passes through the material in the 
jet above and below the accretion disk when the secondary 
passes in front of the primary at superior conjunction. We 
will refer to the size of the jet by the half opening angle, a, 
or the angle between the wall of the conical jet and the jet's 
axis, as shown in Fi gure [T] 

As discussed in iThomas etal] l|201ll 'l. emission features 
can originate anywhere in the RR, whereas absorption fea- 
tures must arise from material in the line of sight to the 
primary star. Therefore, when studying absorption features 
one must use i^s = 35°. Furthermore, the absorbing mate- 
rial most likely lies between the star and the inner edge of 
the circumbinary disk. Our focus for this study will be on 
the periodically varia ble blue- shifted Ha absorption feature 
attributed to the jet l|Witt et al.„200a '). 



2.2 Orbit 

The orbit of the bin ary has been stud i ed pr eviously by 
IWaelkens et al. I l|l996t ) and lThomas et all (|201ll '). the orbital 
period was found to be ~ 317 days. The system is receding 
from Earth at roughly 19 km s~^. In all of our data the obser- 
vations are phase folded using the orbital period of 317 days 
with a zero phase point to = JD 2448300, which was cho- 
sen to correspond roughly to minimum in the radial velocity 
curve IWaelkens et al.ll 19961 ). In this paper all orbital phases 
indicate the location of the primary. The location of the sec- 
ondary is always opposite relative to the centre of mass. To 
aid in the interpretation of the data and model calculations, 
Figure [2] shows the orbits of both stars. We will frequently 
refer to the following key orbital locations: inferior conjunc- 
tion (IC) when the primary is between the observer and the 
secondary, which occurs at phase 4> = 0.21; superior conjunc- 
tion (SC) when the secondary is between the observer and 
the primary (<j!> — 0.59); apastron (0 — 0.79); and periastron 
{<f> = 0.29). 
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Figure 1. Geometrical model of the RR. Panel (a) illustrates 
the overall geometry of the RR as seen from the plane of the 
sky. The height of the circumbinary disk and the inner radius are 
to-scale. The binary at the centre appears 10 times larger than 
the scale of panel (a). The outer edge of the circumbinary disk 
(1850 au) is not shown. The cone walls indicated are meant to 
correspond to the X-shaped structure seen in the RR. The polar 
axis shown is perpendicular to the orbital plane, which is tilted 
with respect to the direction of the observer by 5° (i = 85°). 
The outer nebula is seen via the direct line of sight, while the 
binary is seen under at the effective inclination, i^fs = 35°, via 
the indirect line of sight. Panel (b) is an enlargement of central 
portion of panel (a). Panel (b) shows the indirect line of sight 
passing through the jet at the phase of superior conjunction (SC), 
= 0.59 (see Figure [2J. The size of the primary star and the 
separation between the primary and secondary are shown to-scale 
with 1 au as indicated. The accretion disk is shown around the 
secondary, with the jet emanating from above and below the hot 
accretion disk. The fiducial mark between the stars corresponds 
to the centre of mass (COM). 
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Figure 2. Orbital plane diagram. The orbit of the primary is 
shown by the solid ellipse, while the dashed ellipse shows the 
orbit of the secondary. The orbital phases are indicated on the 
orbit of the primary. The primary star moves in the direction of 
increasing phase nu mbers, anticlockwise. P hase zero was chosen 
to be consistent with lWaelkens et al.l lll996h : it corresponds to JD 
2448300. The locations of the stars for each of our observations 
are shown by the circles on the orbits. The secondary is always 
opposite the centre of mass (COM), which is located at (0,0). 
The locations of apastron and periastron are at the intersection 
of the orbit with the antipodal line. Inferior conjunction (IC), 
phase (/) = 0.21, occurs where the direction vector to Earth (5° 
below the plane of the page) intersects the primary's orbit on 
the Earth side of the COM. The phase of superior conjunction 
(SC), phase (/> = 0.59, occurs where the Earth vector intersects 
the primary's orbit on the far side of the COM. 



3 OBSERVATIONS 

The high velocity biconical jet is observed via absorption 
in the Ha hne profi les. Since the original study of the jet 
bv lWitt et all l|2009l ) 16 new spectra have been obtained in 
order to study the geometry of the jet, with emphasis on 
the orbital phases during which the material in the jet lies 
in the line of sight to the primary star. The data set now 
contains a total of 33 observations, the locations of which 
are shown in Figure [S] 



3.1 Details of the Observations 

The spectra were acquired at the Apache Point Observatory 
(APO ) with the ARCES echelle spectrograph (|Wang et all 
l2003l) and the 3.5-m telescope. The spectrograph has a re- 
solving power of R = 38,000, producing a velocity resolution 
of 8 km s~^. The signal-to- noise ratio of the spectra in the 
Ha region is roughly 1000. Each observation is the average of 
two spectra acquired consecutively on the same night. The 
dates of observation were chosen for the purpose of study- 
ing the periodic variability of t he RR. The data redu ction 
has remained unchanged since iThorburn erabl l|2003l ) and 
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iHobbs et al] ||2004|). a nd the complete data set is described 
fully in lThomas et all l|201lh . 

The Ha line profiles, shown in Figure|31 have several dis- 
tinct features. First, there is a narrow, symmetric, stationary 
emission peak that has a FWHM ~20 k m s~^ l|jura et alj 
ll997l : lHobbs et al.ll2004l : IWitt et al.ll2009D . Second, there is a 
(lower amplitude) broad emission plateau (FWHM of nearly 
200 km s~^). Third, there is an asymmetric blue-shifted ab- 
sorption feature on the blue side of the emission plateau that 
appears shortly after the phase of IC and becomes strongest 
around the time of SC, see panel b of Figure |3] There are 
also broad emission wings that extend to ± 600 km s~^. 
The narrow emission feat ure has been ascribed to material 
in the central Hii re gion l|jura et al.ll 19971 : 1 Witt et al.ll2009l : 
omas et al.ll201lb . The likely source of ionising photons 
is the UV prod uced in the inner regions of the secondary's 
accretion disk (|Witt et al.| [2009). The width of the narrow 
emission feature is similar to the orbital velocities for the in- 
nermost circumbinary material, which suggests that this ma- 
terial is in circumbinary orbit. The broad emission plateau 
and extended emission wings are likely due to either the 
material in the jet, and/or material in the accretion disk. 
Since emission can originate from anywhere in the system 
it is harder to pinpoint its origin. The absorption feature 
has been attribu ted to absorption by material in the jet 
l|Witt et al.ll2009l ). It is this absorption feature we shall fo- 
cus on. 



3.2 Photospheric subtraction 

Analysing the Ha profiles in detail to study the jet 
requires that the observed spectral profiles have the 
contribution of the primary's photosphere removed. 
The NLTE stellar atmosph er e cod e , phoenix, version 
15.04 .0 0E dH auschildt' 1992', 'l993'; Hauschildt fc BaronI 
19951: lAllard fc Ha uschildt 1995.: Baron et al.l 1 19961: 



Hauschildt et al.l 1 19961: 
19971: iBaron fc Hauschildtl 



Hauschildt. Baron, fc AllardI 
19981 : lAllard et al.l l200ll : 



Hauschildt. Lowent hal. fc BaronI 1200 ll ). was used to 



calculate the expected Ha absorption line profile of the 
primary star. Other model lines produce d by phoenix 
are m uch narrower than the observed lines (|Thomas et al.l 
l201lb : this feature likely applies to the Ha line as well. 
Therefore, the subtraction of the phoenix model line from 
the observed line profiles was carried out us i ng th e sub- 
traction technique discussed in iThomas" et all (|201ll ). This 
technique removes the underlying photospheric absorption 
line by accounting for the intrinsic broadening that seems 
to be present in all unblended atomic absorption lines in th e 
spectra of t he RR llHobbs et al.l[20o3: iThomas et aLllioTH l. 
The work of ICrinin. Mitskevich. fc Tambovtseval ()2006l l and 
iGrinin. Tambovtseva. fc Weigeltl l|2012h attributes similar 
broadening to light scattering off moving grains. The model 
line profile is broadened by a broadening kernel derived 
from the observed characteristic photospheric absorption 
profiles. Examples of the subtracted spectra are shown 
in Figure [S] these same spectr a are pres e nted in their 
unsubtracted form in Figure 5 of I Witt et all l|2009l ). 



4 MEASUREMENTS 

In order to study the jet geometry, the equivalent width 
of the blue wing of the broad emission plateau of the Ha 
profile, VFA(Ha) — Wua, was measured as a function of or- 
bital phase. The orbital phase dependence of these mea- 
surements helps to constrain a, the half opening angle of 
the jet. The emission equivalent width was measured for 
each photosphere-subtracted spectrum, from —27 km s~^ to 
—600 km s~^, so as to not include the effects of the cen- 
tral narrow Ha emission feature. The results are shown in 
Figure |4l Emission equivalent widths are typically defined 
to be negative, in Figure |4] we see that the maximum Ha 
equivalent width measured is roughly —3 A. As a function 
of phase we see that the amount of emission decreases just 
after IC, reaching a minimum near SC after which it in- 
creases again as the primary star again approaches IC. The 
observed decrease corresponds to more absorbing material in 
the indirect line of sight, thus lowering the amount of emis- 
sion observed. Since the equivalent width measurements re- 
main negative, the amount of absorption never exceeds the 
amount of emission. The scatter in the data is due to uncer- 
tainties in the continuum level and residual telluric features. 
The equivalent widths were measured by hand and with a 
script, each data point is the average of these two values and 
the error bars represent one standard deviation. 

Another way to visualise the characteristic absorption 
feature of the jet is to construct a difference spectrum by 
subtracting the spectrum that shows no sign of the absorb- 
ing material. The photosphere-subtracted spectrum which 
occurs at IC, phase = 0.21, was taken to be the template 
spectrum due to its symmetric observed profile (see panel a 
of Figure O), and because this is the phase during which the 
emission equivalent width is least affected by the presence 
of the jet (i.e., the equivalent width is most negative), see 
Figure 13] The template spectrum was subtracted from every 
photosphere-subtracted spectrum; the results are presented 
as a dynamical spectrum for all 33 observations in Figure (5) 
The systematic variation in the strength of the absorption as 
a function of orbital phase naturally agrees with the equiv- 
alent width measurements. The majority of the absorption 
occurs within roughly —200 km s~^, while the underlying 
emission extends to approximately —600 km s~^. The ve- 
locity at which the absorption maximum occurs, Vkbs max, 
is relatively constant, roughly —35 km s^^, while the veloc- 
ity of the blue edge of the absorption feature, Vbiuo edge, is 
variable and has a maximum value of roughly —200 km s~^. 



5 MODELING THE JET 

In this section our goal is to determine the jet half opening 
angle, a (the angle between the conical jet axis and the wall 
of the conical jet, see Figure [TJ. To do so we will employ two 
observational diagnostics that depend on the jet geometry. 
The first is the absorption equivalent width, which measures 
the amount of absorbing material in the jet. The second 
is the dynamical absorption spectrum, which measures the 
velocities in the jet. 
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Figure 3. Representative Ho spectra. The logio of tiie flux is plotted on the vertical axis, as a result the continuum level is at logio(l) = 
0. This scaling was done to deemphasize the height of the narrow emission peak. The continuum-normalised spectra have been corrected 
for the underlying photospheric absorption and its intrinsic broadening. Panel (a) shows the most symmetric spectrum, which occurs at 
IC. Panel (b) shows a spectrum at phase <ji = 0.58, close to SC, with its strong blue-shifted absorption. For the relative position of the 
stars at these phases see Figure[2l The lower horizontal axis indicates the centre of mass velocity in km s~^, to convert to heliocentric 
velocity add 19 km s~^. 
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Figure 4. Equivalent width measurements of the blue wing of 
the Ho profile. Larger negative values indicate stronger emission, 
and correspondingly less absorption by the jet. The error bars 
are one standard deviation for the average of two measurement 
techniques. The vertical lines correspond to the phases of IC and 
SC, see Figure[2l The absorption maximum of the jet occurs near 
SC, while the minimum of the absorption occurs near IC. 



Figure 5. Dynamical spectrum of the Ho absorption feature. 
Each spectrum shown is the difference (observed — template). 
The continuum corresponds to 0.00 on the linear grey-scale bar, 
absorption corresponds to negative numbers. The vertical lines 
correspond to the FWHM of the narrow emission peak, any sub- 
traction here is inaccurate due to variations in the intensity of 
the narrow emission peak. 
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5.1 Modeling the equivalent width measurements 

From the observed phase coverage {(p = 0.4 to <j> = 0.9), 
where the primary star is occulted to varying degrees by the 
jet (Figure [4] and Figure [5]), we can say that a must be large 
enough that the jet attenuates light from the primary star 
for most of the orbit. This implies that a is near but does 
not exceed the effective inclination of the line of sight to the 
primary star. If a were larger than ieff there would always 
be some occultation of the primary star. 

It is useful to look at the variation in the equivalent 
width of the absorption feature, instead of the emission fea- 
ture. We used the change in equivalent width with respect 
to the symmetric spectrum at IC. The change in equivalent 
width produced by the absorption feature was calculated by 
subtracting the equivalent width at IC from each equiva- 
lent width measurement (Wifa(<^)— W^). The results were 
normalised for comparison to the normalised model curves 
(to be discussed below) and are shown in Figure |6l In the 
figure we see that as a function of phase the line of sight 
starts to pass through the jet just after IC. The line of sight 
reaches absorption maximum near SC, and trails off toward 
zero absorption as the primary star again approaches IC. 

The path length through the jet was used to approx- 
imate the equivalent width, W\, which on the linear por- 
tion of the curve of growth is proportional to the column 
density, A'^. At constant density, A'^ is proportional to the 
geometric path length, I, through the jet. This simplistic 
model allows us to calculate the orbital phase dependence 
of the normalised path length through the jet at various 
values of a, for comparison to the equivalent width mea- 
surements. Model calculations for three values of a (33°, 
35°, and 37°) are over-plotted on the normalised equivalent 
width measurements in Figure (6] The selected curves are 
calculated for an indirect line of sight angle of icft = 35°, 
the "effective" inclination angle. The model curves fit sur- 
prisingly well for such a simple model. The results indicate 
that the best fitting model is such that a ~ 35° — ioff. Since 
the value of ics is uncertain, the path length curves were 
calculated for effective inclination angles of 30°, 40°, and 
50°. In all cases, the results were such that the value of a 
required to match the measurements was, respectively; 30° , 
40°, and 50°. Therefore, we can state that a must approx- 
imately equal the effective inclination angle for the system, 
regardless of the assumed effective inclination. 

IWitt et al.l l|2009l ) noted that the equivalent width mea- 
surements are asymmetric about the phase of SC. This 
asymmetry was interpreted as a sign of material trailing the 
jet. However, we find that trailing material is not required to 
explain the observed asymmetry in the orbital phase depen- 
dence of the absorption feature. The curves are asymmetric 
due t o the eccentricity of the orbit (e ^ 0.3) (|Thomas et al.l 
[ioil). The eccentricity of the orbit causes the indirect lines 
of sight to the primary star to pass through different heights 
above the orbital plane as a function of phase. Using our 3D 
model we calculate the height above the orbital plane for 
lines of sight to the upper and lower limbs of the star, as 
well as the centre of the primary star (see Figure O. The 
large size of the primary star spreads out the curves for 
lines of sight to the leading and trailing limbs of the star as 
well; however, the change in phase coverage is minor (about 
±0.025 in phase), the curves are otherwise similar. The re- 




0.4 0.6 
Phase 

Figure 6. The data points are the normalised difference in the 
equivalent width measurements for the blue-shifted absorption 
feature of the Ha profile, shown in arbitrary units. The curves 
correspond to the normaUsed geometric path length through the 
jet for different values of a, also in arbitrary units. The dashed, 
solid, and dash-dot curves correspond to a = 33°, 35°, and 37° 
respectively. The horizontal axis corresponds to the orbital phase. 
The vertical lines correspond to IC and SC, see Figure [2] 



suiting path lengths {I — exit point (45 au) — entry point) 
are also asymmetric about SC. Therefore, we need not in- 
voke trailing material to explain the observations. Any such 
trailing material would be of significantly lower density than 
the jet, and would have a significantly smaller effect on the 
observed spectra. If present, this material may simply con- 
tribute to a more or less uniform background present in all 
spectra, and would be effectively removed in the subtraction 
of the template spectrum. 



5.2 Modeling Jet Velocities 

Here we use the dynamical absorption spectrum (Figure O 
to determine the jet kinematics by fitting the observed ve- 
locities. We will focus on the observed velocity of the blue 
edge of the absorption feature, Vbiue edge, and the observed 
velocity of the absorption maximum Vabs max- We will use 
the majuscule, V , to refer to velocities in the data, and the 
minuscule, v, to denote computed velocity values. For our 
calculations below, we used a = iea = 35°. 

We modelled a constant velocity jet, an accelerating jet, 
and a jet with a latitudinal velocity structure. The constant 
velocity model has a velocity of 200 km s~^. The accelerat- 
ing jet had only a radial dependence to the velocity of the 
form v{r) = vo {r /rnorm), where vq is the velocity along the 
jet axis (set to 200 km s~^), r is the radial distance from 
the centre of mass of the secondary star, and r„orm is the 
normalisation radius in this case we used 1 au. In all the 
models we assume the outfiow is purely in the radial direc- 
tion, centred on the secondary star. 
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Figure 7. Height above the orbital plane at which the indirect 
line of sight passes through the jet as a function of phase. Three 
lines of sight are shown. The black curves are for the line of sight 
to the centre of the star. The upper and lower grey curves corre- 
spond, respectively, to the upper and lower limbs of the star. The 
long dashed curves correspond to where the lines of sight enter 
the jet; the solid curves correspond to their intersection with the 
plane in the plane containing the axis of the jet. The vertical lines 
correspond to IC and SC, see Figure [2l 



The best fitting model was the latitudinal velocity 
structure, which has a velocity gradient (with the higher ve- 
locities along the jet axis) in the latitudinal (polar) direction 
only. Typically, acceleration of the outfl owing jet normally 
happ ens within 1 au for solar type stars (|Appenzeller et ahl 
l2005h . after which the motion is ballistic. Since our line of 
sight never probes lower than roughly 1 au, see Figure [71 we 
will ignore such acceleration in this model. Typically, the ter- 
minal speed along a given streamline scales with the escape 
speed at the base of the streamline where the material is 
launched. For this reason hydrodynamic models predict ve- 
locity distributions with the highest velocity material along 
the cone axis (sometimes referred to as a spine or jet) and 
the slower material at t he edges (sometimes ca lled the sheath 
or conical wind) fe.g.. iRomanova et al.ll2009l '). The velocity 
of the launched material is dependent upon distance from 
the centre of mass of the secondary star. We model this de- 
pendency as a latitudinally dependent power law velocity 
structure. 



v{e)f ■ 



Vo + {Va — Vo) 



0' 



(1) 



where 9 is the polar latitudinal coordinate (which varies 
from 0° to a), vq is the velocity along the jet axis (set to 
200 km s~^), Va is the velocity at the edge of the jet (set to 
30 km s~^). By symmetry v{6) must be an even function of 
0, so we set p to an even integer. A constant velocity jet is the 
case p = 0, and for the spine sheath case we set p = 2. Figure 
[8]illustrates a simplified system geometry indicating the unit 
vectors and coordinates used in this model. As a function 



Indirect 
line of sight 




Axis 



Figure 8. Simplified jet diagram illustrating the velocity struc- 
ture described in the text. The sketch is not to-scale. For sim- 
plicity, only one indirect line of sight is shown, the unit vector 
n describes the direction of the line of sight. The latitudinal po- 
lar coordinate, 9, is used to parametrise the velocity structure in 
Equation [T] All velocities are directed radially, f. The entry and 
exit points occur at the locations indicated at the intersections of 
line of sight with the jet wall. 



of orbital phase, projected line of sight velocities along the 
vof ■ h and jet entrance v]f'^ = Vaf ■ h are 



jet axis Uo°^ 



over-plotted on the dynamical spectrum in Figure [9] The 
maximum of uo°^ is roughly 160 km s~^ and the orbital 
phase dependence of the projected jet axis velocity matches 
the blue edge of the absorption feature, Vbiuo edge, reason- 
ably well. The calculations show that projected jet entrance 
velocity w^^^ is roughly constant at 35 km s~^, which is 
where the absorption maximum T4bs max occurs in the data. 
We find that this simple model fits the data surprisingly well, 
matching not only the observed velocities, but the observed 
orbital phase dependence as well. Any arbitrary model can 
fit a single observation, but the fact that one model fits all 
of the observed data so well makes this fit reasonable. In- 
cluded in Figure |9] are curves representing the other two 
models (grey curves). The grey curves are for the same lo- 
cations in the jet as the best fitting model. It is clear that 
the grey curves do not fit the data. 
It appears that Vabs max and 

are related to each 

other, and that Vbiue edge and 

t;5'OS are also related to one 
another (the best fit to the data gives vo — 200 km s"^ and 
Va — 30 km s~^). To understand this relation, let us look at 
why it might be the case. Continuum photons from the pri- 
mary star can only be absorbed via a transition line (e.g.. 
Ha) if the frequency of the photon lies within the range 
of frequencies spanned by the line profile. However, in the 
jet the material is moving. Therefore, the transition line in 
the jet must be Doppler shifted to the velocity of the ab- 
sorbing particle. The continuum photon from the star can 
only be absorbed if its frequency lies within the line width 
(Avo) frequencies of the Doppler shifted transition line. The 
geometric length, ds, through the jet in which the contin- 
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Figure 9. The dynamical spectrum shown is identical to Fig- 
ure [5] The overlain black curves correspond to best fitting model 
(latitudinally dependant). The solid black curve is the projected 
line of sight velocity (PLOSV) along the jet axis, i'q'"'^, while the 
dotted black line corresponds to the PLOSV at the entrance to 
the jetjiJ^*^^. The grey curve with square markers is the w^'-'^ for 
the constant velocity model, while the grey curve with triangle 
markers is the ■Ug for the accelerating case. The grey curve with 
circle markers is the ii^*-*^ for the constant velocity case, while the 
grey curve with the upside-down triangle markers is the v^*-"^ for 
the accelerating model. The range of phases from <f> = 0.05 to 
(f) = 0.3 show no curves because the indirect line of sight does not 
pass through the jet. 

uum photon may be absorbed depends on the width of the 
line and the velocity gradient (dvLOS /ds) of the moving ab- 
sorbers. In the limit that the line is very narrow, the ge- 
ometric path length (Sobolev length) depends only on the 
velocity gradient. A large gradient in the velocity of the 
absorbing material in the jet will lead to a narrow interac- 
tion region, while a smaller velocity gradient gives rise to a 
longer interaction region. The longer the interaction region 
the stronger the absorption. Therefore, we expect that the 
absorption will be strongest near the walls of the jet, due to 
the relatively small velocity gradient {Av{6)f / A9) and cor- 
respondingly longer interaction length. We expect that the 
material near the jet axis will have smaller absorption due to 
the larger velocity gradient, while the material near the en- 
try point will have intermediate absorption values. We now 
have a physical explanation of why Vkbs max = v]^'^^ , and 

1/ — ,,LOS 

vblue e dge — t^Q 

In IWitt et ai] (|2009l 'l it was noted that the velocity of 
the outflow was variable based on the observed variation 
of Vbiuc edge , see Figure [9] This variability was interpreted 
as a jet that varies in launch speed or size as a function 
of orbital phase. The range in the variation of the veloc- 
ity would imply that the disk is substantially changing size, 
since the velocity of the launched material depends on the 
distance from the centre of the secondary star. Such a dra- 
matic change in the radius of accretion disk is in contradic- 



tion to the fall-in time as calculated in IWitt et al.l (|2009l ) of 
roughly 500 years. The fact that the accretion disk is stable 
for longer periods of times than the orbital period implies 
that the jet is likely stable as well. Variable accretion rates 
could possibly affect material launched from the edges of 
the disk; however, this material should be launched with 
slower velocities and cannot explain the observed variation 
in Vbiue edge- In contrast, as shown by our calculations, a 
conical jet with a latitudinally-dependent velocity structure 
and fixed half opening angle angle, a, fits the observations 
remarkably well, and our model does not require a complex 
interpretation involving a jet that changes size or launch 
speed with orbital phase. Finally, we note that neither a 
constant velocity jet (p = 0) nor a radially accelerating jet 
(without angular variation) can reproduce the phase depen- 
dence of the velocities of the blue edge, Vbiue edge, and the 
absorption maximum, Vkbs max. 



6 DISCUSSION AND CONCLUSIONS 

In this study we analysed the blue shifted absorption fea- 
ture previously attributed to the bipolar jet emanati ng from 
a hot accretion disk around the secondary star ( Wi tt et alj 
'2009!). This jet system shares many similari ties with jet ob- 
served in BD-|-46°442 (jCorlova et al.ll2012l ). Both systems 
are evolved binaries with circumbinary discs, and both sys- 
tems have wide jets, which are launched from the secondary 
stars. However, the viewing geometry of the two objects dif- 
fer. 

Based on the observed range in orbital phase during 
which the primary star is occulted, we find the half opening 
angle of the jet, a, to be slightly less than the effective incli- 
nation angle (35°). It cannot be greater than this angle, as 
the spectrum is virtually unaffected near the phase of infe- 
rior conjunction (IC). This conclusion is ind e pende nt of the 
assumed inclination angle. IWaelkens^ et al.1 l|l996t ) showed 
that they needed an effective inclination of 35° in order to 
fit the observed periodic variation in the visual magnitude; 
however, they used the geometry of the X-shaped bi-cone 
to infer this effective inclination angle. That said, since our 
jet half opening angle is essentially the same as the viewing 
angle it seems very likely that the observed jet is actively 
carving out the biconical cavity in the otherwise spherical 
outfiow from the AGB stage of the primary star's life. Fur- 
thermore, in order to fit the equivalent width measurements, 
a wide jet is required. Therefore, our observations are in- 
consistent w ith a narrow precessin g jet model, as recently 
proposed by I Velazquez et al] l|201ll ). 

The velocity of the blue edge of the absorption, 
Vbiue edge , IS Seen to Vary. We conclude that this variability is 
likely the result of our indirect line of sight passing through 
a jet with a latitudinal velocity structure, and not to a jet 
that is variable with orbital phase. The jet is stable over 
the duration of the binary orbit, and does not show varia- 
tion of outfiow velocity or direction. Over the 8 year span of 
observations no precession in the jet is apparent. The max- 
imum outflow velocity in the jet is wo = 200 km s~^ on the 
jet axis, and Vbiue edge is physically related to this velocity. 
The velocity of the absorption maximum, Vabs max, in the 
jet occurs at approximately —35 km s"'^, which is physically 
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related to the velocity at the walls of the jet, which we find 
to be Va = 30 km s~^. 
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